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We study the effects of electron-phonon interactions on the transport properties of a molecular 
quantum dot coupled to two Luttinger-liquid leads. In particular, we investigate the effects on the 
steady state current and DC noise characteristics. We consider both equilibrated and unequilibrated 
on-dot phonons. The density matrix formalism is applied in the high temperature approximation 
and the resulting semi-classical rate equation is numerically solved for various strengths of electron- 
electron interactions in the leads and electron-phonon coupling. The current and the noise are in 
general smeared out and suppressed due to intralead electron interaction. On the other hand, the 
Fano factor, which measures the noise normalized by the current, is more enhanced as the intralead 
interaction becomes stronger. As the electron-phonon coupling becomes greater than order one, the 
Fano factor exhibits super-Poissonian behaviour. 
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I. INTRODUCTION 

In the relentless search for smaller electronic de- 
vices, the idea of fabricating extremely small transis- 
tors using quantum dots has become an important topic. 
In recent years, the possibility of using a very small 
molecule as the quantum dot in a transistor has become 
apparenliiiSi 3 ^^^^. Such a device may be thought of 
as a molecular quantum dot weakly coupled to macro- 
scopic charge reservoirs, or leads. 

Molecules integrated into these transistors often have 
very complex structures and therefore introduce new 
electronic transport properties. For example, the shape 
of the molecule or the position of the molecule relative 
to the leads can change as charges are added onto the 
molecule. A number of experiments have investigated 
the effects of coupling between tunneling electrons and 
the on-dot vibrational quanta on various transport prop- 
erties of the molecule2ii2iiiii2. Of particular interest is 
the experimental evidence for electron-phonon coupling 
in a molecular device composed of Cgo molecule deposited 
between a pair of gold electrodes^. Experimental results 
show peaks in the differential conductance which may be 
due to the effect of coupling between tunneling electrons 
and the molecular vibrational mode. 

Several theoretical worka 13 i 14 i 15 i 16 i 17 i 18 i 19 i 20 have har- 
nessed the study of electron-phonon coupling by intro- 
ducing an on-dot phonon degree of freedom that couples 
to on-dot electrons and observed the effects the coupling 
has on various transport properties of the molecular tran- 
sistor. One of the most recent comprehensive work on 
phonon effects in molecular transistors coupled to two 
non-interacting leads was conducted by Mitra et al*&. 
They developed the density matrix formalism to study 
transport properties for both in-equilibrium and out-of- 
equilibrium phonon distributions. 

In our work, we investigate how the phonon distribu- 
tion and electron-electron interaction in the leads affect 
the steady state DC current and the DC noise character- 



istics of a molecular transistor. The phonon distribution 
can drastically change depending on how they are al- 
lowed to relax as they couple onto a bath. As discussed 
in Mitra et ali&, we must consider two competing time 
scales: time scale for the phonons to reach equilibrium 
via their interaction with the bath, r p h; and the dwell 
time of the electron on the dot, Tdwcii- The two limit- 
ing cases considered in this work are: the phonons are 
equilibrated to the bath corresponding to Td wc ii S> r p h; 
and the opposite case where the tunneling electrons see 
unequilibrated phonons corresponding to r p h ^> Tdweii- 

We are particularly interested in the effects of electron- 
electron interaction in the leads. Up to this date, theo- 
retical works have only considered non-interacting elec- 
trons in the leads. This is sufficient if the leads are two 
or three dimensional electron gases, where interactions 
affect the low energy properties of the lead only pertur- 
batively. However, in one dimension, arbitrarily weak 
electron interaction completely changes the ground state. 
In our work, we consider the leads to be a Luttinger liq- 
uid. Our work may represent a physical device composed 
of a molecular quantum dot, such as Ceo, coupled to two 
metallic single- walled carbon nanotubes, which behave 
essentially as interacting one-dimensional objects. 

The outline of this paper is as follows. We consider 
a molecule with a single level coupled to two Luttinger- 
liquid leads. We suppose that the electrons are coupled 
to the on-dot vibrational mode and consider both equili- 
brated and unequilibrated phonons. We solve the model 
in the density matrix formalism and apply the high- 
temperature approximation. The approximation will be 
used to derive the rate equation for the dot occupation 
probabilities and use the equation to compute the steady 
state DC current and the DC noise characteristics of the 
molecular transistor as a function of source-drain volt- 
age. In section [H] we introduce our model, provide a 
brief outline of the density matrix formalism, and present 
expressions for the tunneling current and the DC noise. 
In section ITTT1 we present numerical results for the cur- 
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rent and the noise and discuss the dependence of these 
transport properties on the phonon distribution and the 
intralead electron interaction. We conclude in section llVl 



II. FORMULATION OF THE PROBLEM 
A. Model 

Our model considers a molecule with a single non- 
degenerate energy level, of energy eo, coupled to two 
Luttinger- liquid leads labelled L (left) and R (right). We 
allow the tunneling electrons to couple to a single phonon 
mode, with frequency loq corresponding to the internal vi- 
brational mode of the molecule. The full Hamiltonian is 
given by 



H = H, 



dot 



Hleads + Ht, (1) 

where Hdot describes the on-dot electrons coupled to the 
single phonon mode, Hi ea ds represents the electrons in 
the leads, and H t corresponds to the tunneling between 
the dot and the leads. 

More specifically, Hdot has the form, 

Hdot = eo Yl c ^ c ° + jYl cX ° c ° ( c « Ca ~ 1 ) 

a U \ a / 

+ Xlu (& f + b) C W + "otfb, (2) 

a 

where c\ is the creation operator of electrons with spin 
a on the dot and b^ is the creation operator of the on- 
dot phonons. A is the electron-phonon coupling strength 
and U is the charging energy of the molecule. Our work 
focuses on the limit U — > oo, which may be relevant 
when considering the case of Ceo transistors whose charge 
states are most likely zero or on o 12 i 26 . Hi ea d s and H t are 
given by 



and 



Hl ea ds — J]] Hi ea d-i — J]] Hkin- 
i=L,R i 



Hint- 



h.c. 



(3) 



(4) 



where 

Hkin-i — v F ^2 _ k F )a\ jj i a kt(T . i 

-{k + k F )b{ „M,o-,-, 



(5) 



Hint—i 

l 
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k,k' , f7.tr' 

+ b l-a,a, l b k^ b l+ q ,a'^ b k'^a] } • (6) 



Hkin-i and Hi nt -i respectively correspond to the kinetic 
and the interaction terms of the electrons in lead i; L s is 
the linear size of the leads, and v F (k F ) is the Fermi veloc- 
ity (wavevector) . Here a\ ai (b\ ai ) creates a right(left)- 
moving electron with momentum k and spin a on lead i. 
<72 and g± represent forward scatterings; in our work, we 
will not consider the back-scattering interaction. 

In the absence of back-scattering, the Luttinger Hamil- 
tonian, Hiead-i, is exactly soluble using the technique of 
bosonizationS&S 4 * 2 ^. First, we rewrite the Hamiltonian 
using the Fourier components of the particle density op- 
erator for right and left movers: 



k 

In the continuum limit (L s — > oo), 

/>oo 

H km ^^v F / dq [p + ,(7,i(g)p+, -,i(~g) 

Jo 



(7) 
(8) 



+P-,<T,i(-9)p-,<r,i(l)] . ( 9 ) 



Hint-i = J dq^2{g 2 p+,a,i{l)P-^',t(- < l) 

f7.tr' 

+y \p+,^i(9)P+,a',i(-q) + p-,a,i(-q)p-,a',i(q)]j ■ 

(10) 

In order to express the Hamiltonian in diagonal form, 
we introduce the canonically conjugate Boson operators, 
6~i and n CT ,i, 



(j><y,i{x) 



1 Jq^O 1 



(11) 



and 



n CT ,i(a;) 
1 

2^ 



q^O 



dqe—r-** [ P+ ,,M ~ P-,*M], (12) 



where a is a small convergence factor or the inverse of a 
large ultraviolet cutoff. These operators are transformed 
into the charge and spin boson operators, 



9 c(s),i 



n 



c(s) ,2 



n T ,, + (-)n i;t 

V2 



which obey 



it i(x), U„,i(y)] = i8^ v 8(x-y). 



(13) 
(14) 

(15) 
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The single fermion operators for right- and left-moving 
electrons with spin a on lead i can be written in position 
representation as 

i>±,<r,i(%) = 

lim 1 e ±itj.i->jg[±(^, i +^,,i)+(fl e , ( +rf.,i)l 
where 



(17) 



The Hamiltonian, Hi ea d-i, can be expressed in terms of 
the canonically conjugate charge and spin boson opera- 
tors fEas ll3H4|) as, 



H, 



lead—i 



where 



?; F + — - — 

7T / V 7T 



7TOf + g 4 - g 2 
7TU F + 5(4 + g 2 

v s = wf; 
= l. 



(18) 



(19) 

(20) 

(21) 
(22) 



Notice that K c < 1 (K c > 1) for repulsive (attractive) 
interaction. 

Now a convenient canonical transformations* can be 
applied on the full Hamiltonian making it diagonal in 
the dot variables. The desired transformation is H — > 
H' = e s He~ s where 



,t 



Then, the transformed Hamiltonian is given by 
H' = H' dot + Hi eads + H t , 

where 

H'dot = E 4* + yE 4* (E 4c - I' 

a a \ cr y 



(23) 



(24) 



wo^fo', (25) 



and 



i=L,R 



The transformed phonon operator is 

&' = 6-aE4c 



(26) 



(27) 



which shows that ground state energy of the phonon spec- 
trum depends on electron occupancy of the dot. The 
shifted energy level of the molecule is e' = e — \ 2 uj , 
while the shifted charging energy, U', is not explicitly 
shown since we only consider the limit U —> oo. 



B. Formalism 

1. The high-temperature approximation 

The spirit of the high-temperature approximation 
stems from the presence of two competing energy scales: 
the temperature, /cbT, and the tunneling rate of an elec- 
tron from a lead onto the dot (or vice versa), hT. In the 
approximation, one assumes that ksT hT. Low tun- 
neling rate infers that the number of charges on the dot is 
a well-defined integer, and thus one can characterize the 
state of the dot by a set of dot occupational probabilities. 
The equation of motion of these dot occupational prob- 
abilities is the rate equation. In fact, a dot-lead system 
can be tuned into the low tunneling regime in certain 
experiment: 
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2. The rate equation 

The central element in determining the transport prop- 
erties of our molecular transistor in the high-temperature 
approximation is the semi-classical rate equation, which 
are the equations of motion for the various occupational 
probabilities of the dot. To obtain the rate equation, 
we use the density matrix formalism, developed by Mi- 
tra et ali&, which allows us to determine probabilities 
for various states of the dot under both equilibrium and 
out-of-equilibrium conditions. 

We start with the full density matrix p that obeys the 
equation of motion 



P = -i[H,p]. 



(28) 



The initial step is to assume that the leads are in equi- 
librium independent of the molecular state. To harness 
this assumption, we express p as the sum of the projected 
density matrix, 



Ps = Tr /eads {p(i)} (g) Pleads 

and the complementary density matrix 
Pt = P- Ps- 



(29) 



(30) 



Here, pi ea ds is the density matrix of the two leads at ther- 
mal equilibrium with chemical potential pj, for the left 
lead and pn for the right lead. Tri ea d s denotes tracing 
over the leads degree of freedom. On physical grounds, 
this decomposition is useful because the diagonal compo- 
nents of the reduced density matrix, pr> = Tri e ads{Ps}, 
relate directly to the occupational probabilities of vari- 
ous states of the dot. In the high-temperature approxi- 
mation, the equation of motion for p s , using Eas l28l2"9l 
can be written as 



Psl{t) ~ -iT^leads[H t l(t), Ptl{t)} <8> Pleads 



(31) 



where subscript I indicates that the operators are in 
the interaction picture: Oi(t) = e lHot O{t)e~ lHot , where 
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Hq = Ho + Hi ea ds- Eq|^ indicates that the time evolu- 
tion of p s is coupled to p t and vice versa. However, one 
can show that the equations of motion for p s and p t can 
be decoupled in the high temperature approximation 13 , 
leading to 

dp 8 {t) _ dip D <g> pi eads ) 



dt 

1 

2 



dt 



dr 



'H u e^ H ^- T \H u psit)Y Ha{t ~ T) ] ■ (32) 



By tracing out the leads degree of freedom on both sides 
of Eq|221 one obtains coupled equations of motion for the 
various occupational probabilities of the dot. Identifying 
P™ as the probability for the dot to be in a state with n 
electrons and q phonons, 



p q = ( n ,q\PD\n,q), 
we arrive at the rate equation, 

P» = ^((q-p^o + Uin-l))^^ 1 

i = L,R 
P 

+ n l ((p-q)cj + Un) T^P^ 1 

-M(g-p)wo + tf(n-l))lt,f? 
-2 ni ((p-q)uo + Un)Tl q P?, 



(33) 



(34) 



where 



Wo 



^ U(q- l)\(l + \p- q \)\ 



(35) 



ni(cj) = 



i r°° 

— / dte 

No J-oo 






e 


IT 








r(f + i< 



2ttT 



0-1 



(36) 



m(cj) = 



w / die*-* (WWV>L(0)> 



2 



2tt 

rf 



(37) 



2ttT 



res) 

where ATg = (u c wf) -1 ^ 2 - r* p denotes the rate at which 
an electron on the dot hops off onto lead i or an elec- 
tron on lead i hops onto the dot while the dot phonon 
occupancy changes from q to p during the process; it 
is symmetric under the interchange of q and p. T(z) is 



the gamma function, while n(oj) is the electron occupa- 
tion number for interacting electrons and is analogous to 
the Fermi distribution, f(u>) for non-interacting electrons 
{n(u!) is analogous to (1 — f(ui))). The factors of 2 in the 
first and the last terms in Eql^Jare due to electron's spin 
degree of freedom. The exponent [3 is given by 



1 



(38) 



In the non-interacting limit, K c = 1 and hence (9 = 1. 
However, for any non-vanishing interaction, K c ^ 1 and 
hence (3 > 1. 

The general form of the rate equation is the same as 
the original rate equation for the case of non-interacting 
leads 13 . However, notice that when the interaction in the 
leads is taken into account, the Fermi functions, /(^), in 
the original rate equation, is replaced by the Luttinger- 
liquid distribution function, n(uj). 

In the limit U — > oo, n in Eq|M|is either or 1. Dur- 
ing numerical calculations, we allow q to take on values 
between and some large cut-off, q ma x ^> 1. In this case, 
P™ is a vector of length 2q max , and the rate equation can 
be expressed in matrix form, 



P = MP, 



(39) 



where M is a 2q max x2q max matrix. At steady state, P = 
0, so the solution for P is the eigenvector corresponding 
to the zero eigenvalue of M.. For equilibrated phonons, 
we use the ansatai^, 



P™ = p n e - quJo/T {l - e^ o/T ). 



(40) 



where P° is the probability that the dot is empty and 
P 1 = 1 - P°. 



3. I- V characteristics 

In the density matrix formalism, current through lead 
i is given by 



(I i ) = Tr{p(t)I i } = Tr{p t (t)I i }, 
where the current operator is 



(41) 



h = it, [e- A(bt - fc) (< fc , CT + h\, k Jc a h.c] . (42) 

k.a 

Performing the trace over the leads and the dot degrees 
of freedom, we can express the steady state current in 
terms of dot occupational probabilities, P™. In the limit 

U — ► oo , we get 

{h) =2[2i*n i ((p-g)wo)r* J , 

<3,P 

-Pfrdq-p)^)^]. (43) 
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4- DC noise characteristics 

It is also interesting to see how the second moment 
in current is affected by the electron-phonon coupling 
and electron-electron interactions in the leads. Current 
through the molecular transistor fluctuates even when 
the system is subject to a DC bias. These current fluc- 
tuations, or noise, can be characterized by the zero fre- 
quency component of the Fourier transform S(u>) of the 
current correlation function, 



S LL (t) = \([SI L (t),SI L (0)} + ) 

= ±Tr{p(I L (t)I L (0) + I L (0)I L (t))} 

-m P i L }} 2 , 



(44) 



where 5Ii(t) — Ii(t) — (Ii(t)). Charge conservation im- 
plies II = —Ir- Therefore, the current fluctuations are 
equal at both leads, giving SllW = Snn(t). Therefore, 
it suffices to compute the current correlation only at the 
left lead. 

In the Heisenberg picture, current evolves over time as 



per I L (t) 
we get 



iHl 



I L (0)e 



iHt 



Substituting this into Eal4*4l 



S LL (t) = -Tr{R{t)I L (0)}, 



where 



R(t) = e~ tllt [6I L (0)p + p6I L (0)] 



jm 



(45) 



(46) 



Since SllW = Sll{— t), Sll^) w iU be computed only 
for positive time. The equation of motion for the causal 
function R(t) can be obtained from the equation of mo- 
tion for the density matrix and is given by 

R(t) = 

-i[H,R(t)]+S(t)[6I L (0)p + pSI L (0)}, t>0, (47) 

and is zero elsewhere. Solving the equation of motion 
for R(t) is similar to the procedure for solving the den- 
sity matrix equation of motion (section III B 2|) . First, we 
decompose the causal function into the diagonal compo- 
nent in the dot and leads variables, and the off-diagonal 
component in those variables: R = R s + R t , where 
R s = Rd® Rieads with R D = Tn eads {R s } . From EqH2 
we obtain the equation of motion for R s (t) in the high- 
temperature approximation, 

R„(t) « -i[H t , R t ] + 5(t) (I Pt + Pt I - 2(1) p,) . (48) 



Notice that the time evolution of R s couples to R t . The 
coupled equations of motion for R s and Rt can be decou- 
pled in the high temperature approximation. Then, as 
was done for the density matrix, we trace out the leads 
degree of freedom in the decoupled equation of motion 
for R s and arrive at the equation of motion for Rd. 

R D (t) = MR D + S{t)h, (49) 
where M. is the same matrix as the one appearing in 
Eg and vector h is given by 

K* = -2(I L )P- 

+ 2J2 [(3 - V((<Z ~ pVo + U(n ~ l))T^ p 

p 

-(n + l)P; +1 n L ((P - q)^o + Un)Tl p ] . (50) 

The final step is to express the DC noise in terms of dot 
occupational probabilities and the components of Ru. 
We first rewrite Eg 1451 as 



S LL (t) = ^Tr{R t (t)I L (0)}. 



(51) 



The trace over the dot and the leads degrees of freedom 
can be performed by following the procedure used for the 
density matrix calculation. The resulting expression can 
be Fourier transformed via, 



S ll {lo = 0) = 2 dtS LL (t) 



(52) 



and we arrive at the zero frequency current noise (U — » 

oo), 

S ll (cj = 0) = 2^2 [2(R(0)° q +P°)n L ((q-p)u; )rL p 
i,p 

-(R(0)l - P^)fi L ((q-p)cj )T^ q ] , (53) 



where R(0) is the Fourier transform of vector R_d(£) at 
lu = and R/j(t) is given by 

R q l (t) = (n,q\R D (t)\n,q). (54) 

Here R(0) can be obtained, at steady state, via Eol49l 

R(0) = -M^h. (55) 
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FIG. 1: Tunneling current plotted as a function of source-drain voltage with A = 1.0. E and NE stand for equilibrated and 
non-equilibrated phonons respectively, and S and A stand for symmetric and asymmetric biases respectively. /3 is the exponent 
in the Luttinger-liquid distribution, which, in the noninteracting limit, is 1 and bigger than 1 for both repulsive and attractive 
interactions. 



III. RESULTS 



A. I-V characteristics 



The tunneling current is obtained at steady state for both equilibrated and non-equilibrated phonons. Initially, 
the dot occupational probability vector, P™, is numerically obtained via the rate equation fEa !39[) . For equilibrated 
phonons, the ansatz, Eg 141)1 is used. The current is then numerically computed using Eg 1431 

The current is plotted for three different electron-phonon coupling constants A = 0.5, 1.0, 2.0. For each A, current is 
computed under symmetric bias (fiL = — fin — V s d/2) and completely asymmetric bias (fin = 0, fiL = V s d) conditions 
and both equilibrated and unequilibrated phonon distributions are considered. Each plot consists of five plot-lines 
with different line styles. Each line style corresponds to a particular value of the exponent, f3, in the Luttinger-liquid 
distribution function. The correspondence is shown in Fig^a). All other similar plots follow the same convention. 
All current plots are produced at T = 0.05o;o- 

We first compare currents for equilibrated and non-equilibrated phonons. Mitra et ali& pointed out a peculiar 
behaviour in the current; near A = 1.0 and with symmetric bias, the current is larger for equilibrated phonons 
than for non-equilibrated phonons. This is an unusual result since one would expect that phonons would arrange 
themselves so as to maximize the current in the non-equilibrium case. They attribute this peculiarity to the choice 
of A at which many of the higher-order diagonal matrix elements (corresponding to g-phonon-g-phonon processes), 
Li , are suppressed. Although their work involved non-interacting leads, one can see in FigUthat this behaviour is 
observed also in the case of interacting leads. The behaviour, however, is not observed when A = 2, consistent with 
Mitra et aVs claim that this anomalous behaviour is pertinent to A values close to 1. 

The most apparent consequence of introducing interaction in the leads is the suppression of the tunneling current; 
the current decreases as interaction strength is raised. The question is whether the suppression of the current is 
mainly due to changes in the phonon distribution induced by electron interaction in the leads or a more direct 
consequence of the Luttinger-liquid distribution function. Fig^ plots the phonon probability distribution for A = 1.0. 
Upon careful inspection, a slight narrowing in the phonon distribution is observed; probabilities of dot states with 
low phonon numbers are increased while the wing of the distribution, corresponding to the dot states with high 
phonon numbers, is suppressed. Although they are not shown, the same subtle change in the phonon distribution 
is observed for other values of the electron-phonon coupling. However, because these changes are very small, the 
suppression in the current is unlikely due to the changes in the phonon distribution. Thus the main effect comes from 
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FIG. 2: Tunneling current plotted as a function of source-drain voltage with A = 0.5. 
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3: Tunneling current plotted as a function of source-drain voltage with A = 2.0. 



the Luttinger-liquid distribution function that has replaced the Fermi distribution in the expression for the current 
fEci 14311 . The Luttinger-liquid distribution function represents single-particle- like excitations decaying faster over time 
than non-interacting electrons. This results in the suppression of the overall current. 

Although Fig0|shows that the phonon distribution is not a strong function of electron interaction in the leads, Mitra 
et alJ2> showed that the phonon distribution is strongly affected by the electron-phonon coupling, A. In particular, 
the distribution was farther from equilibrium for weak coupling than for strong coupling. We have observed that 
the phonon distribution strongly depends on the electron-phonon coupling even in the presence of intralead electron 
interaction. We conclude that varying the electron-phonon coupling has a much larger effect on the phonon distribution 
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FIG. 5: DC noise plotted as a function of source-drain voltage with A = 1.0. 



than varying the interaction in the leads. 

Another noticeable effect on the tunneling current due to interaction in the leads is the smearing and change in 
the slopes of the plateaus. These plateaus, which occur in a range of bias where the tunneling electrons encounter no 
new accessible excited state of the molecule, are flat in the case of non-interacting leads. However, as the interaction 
increases, the slopes also increase. This observation can be explained qualitatively in the zero-temperature limit, which 
is not the appropriate limit for this work but may give insight into this issue. The single-particle time correlation 
function for the Luttinger-liquid behaves as l/t@ . Fourier transform of this correlation function, which relates directly 
to n(u>), then behaves as io^~ x . Therefore, the current should show V^ -1 -type behaviour at low voltages. Each time 
you excite a new phonon, current should show (V s d — Vb)' 3_1 -type behaviour, where Vb is the voltage at which a new 
phonon is excited. Indeed, in the non-interacting limit, where (3 = 1, the slope is zero. However, when /3 > 1, the tilt 
and the smearing of the plateaus result. 



B. DC noise characteristics 



DC noise characteristics are obtained under both symmetric and asymmetric bias conditions for non-equilibrated 
phonons. The noise for equilibrated phonons is not presented here. In section All Al we argued that the modification 
in the current due to interaction in the leads is mainly due to the Luttinger-liquid nature of the electrons in the 
leads, not so much the induced change in the phonon distribution. For this reason, we anticipate that the effects of 
interaction on the DC noise for equilibrated phonons will be very similar for the case of non-equilibrated phonons. 
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FIG. 6: DC noise plotted as a function of source-drain voltage with A = 0.5. 
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FIG. 7: DC noise plotted as a function of source-drain voltage with A = 2.0. 



The noise is plotted at the same three values of A. The correspondence between the plot-line styles and the 
interaction parameters follow the same convention as in the plots for the current (see Fig.l). Figs l5l7l show that the 
noise is suppressed by interaction in the leads. This is expected since interaction in the leads causes tunneling events 
to be more correlated. The increase in correlation among subsequent tunneling events leads to decreased fluctuations 
and hence decreased current noise. 

In Figs 1811(11 we plot the current-normalized DC noise, or the Fano factor F, 

mv , S(V sd ) - S(V sd = 0) 

F{Vsd) = , r • (56) 

2el(Vsd) 

For all values of A, electron interaction in the leads results in an increase in F. The plots also show that F depends 
dramatically on A. In particular, F grows rapidly as A increases^!, and is found to be super-Poissonian for A > 1 (see 
FigHEJ. 

Fie ll II plots the noise at a higher temperature, T = O.lwo- Raising the temperature has a similar qualitative effect 
on the noise as increasing the interaction in the leads; the steps in the noise curve are smeared out and the plateaus 
between steps acquire a slope. However, the overall noise amplitude does not change when the temperature is raised 
while the amplitude decreases when the interaction is introduced. 



IV. CONCLUSION 

In this paper, we investigated the steady state current 
and the DC noise characteristics of a single level quan- 
tum dot coupled to two Luttinger liquid leads. In partic- 
ular, we studied how the electron interaction in the leads 
and the coupling between tunneling electrons and on-dot 
phonons affect the dot's transport properties. We con- 
sidered both equilibrium and out-of-equilibrium phonon 
distributions, and both symmetric and asymmetric bias 



orientations. The density matrix formalism was used to 
harness both phonon distributions in our calculations and 
to derive the semi-classical rate equation for the dot oc- 
cupational probabilities. The calculations were done in 
the high-temperature approximation, in which the elec- 
tron tunneling rate was assumed to be small enough so 
that the dot charge fluctuations are insignificant. 

We have found that interaction in the leads in general 
suppresses the overall tunneling current and the noise, 
and smears out the steps in the I-V characteristics. This 
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FIG. 9: Fano factor plotted as a function of source-drain voltage with A = 0.5. 



effect is mainly a consequence of the Luttinger-liquid cor- 
relation in the leads. The change in the phonon distri- 
bution of the dot due to electron-electron interaction in 
the leads is relatively small albeit we can see a slight 
narrowing of the phonon distribution as the interaction 
strength is increased, leading to only a minor contribu- 
tion to the changes in current and noise. Comparing 
these behaviours of the current and noise with the ef- 
fect of increased temperature, we find that increasing the 
temperature does not suppress overall current and noise 
even though it smears out the steps in I-V and noise 



characteristics. It is also found that the Fano factor, 
that represents the noise normalized by the current, be- 
comes enhanced as the intralead electron interaction gets 
bigger. The Fano factor increases rapidly as the electron- 
phonon interaction becomes stronger and, interestingly, 
shows super-Poissonian behaviour as the electron-phonon 
coupling becomes greater than order one. 
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